Analytical treatment of SUSY Quasi-normal modes in a 
non-rotating Schwarzschild black hole. 

1,2 Pablo Alejandro Sanchez*, 1,2 Jesus Martm Romero*, 1,2 Mauricio Bellini* 

1 Departamento de Fisica, 
Facultad de Ciencias Exactas y Naturales, 
Universidad Nacional de Mar del Plata, 
Funes 3350, (7600) Mar del Plata, Argentina. 

2 Instituto de Investigaciones Fisicas de Mar del Plata (IFIMAR), 
Consejo Nacional de Investigaciones Cientificas y Tecnicas ( CONICET), Argentina. 

Abstract 

We use the Fock-lvanenko formalism to obtain the Dirac equation which describes the interaction 
of a massless 1/2-spin neutral fermion with a gravitational field around a Schwarzschild black hole 
(BH). We obtain approximated analytical solutions for the eigenvalues of the energy (quasi-normal 
frequencies) and their corresponding eigenstates (quasi- normal states). The interesting result is 
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energy eigenvalue. 
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I. INTRODUCTION 



The BH can be understood as a thermodynamical system whose (Hawking) temperature 
and entropy are given in terms of its global characteristics (mass, charge and angular mo- 
mentum). They are obtained by solving a wave equation for small fluctuations subject to 
the conditions that the flux be ingoing at the horizon and outgoing at asymptotic infinity. 
Generally, these conditions lead to a set of discrete complex eigenfrequencies, where the real 
part represents the frequency of oscillation and the imaginary part represents the damping. 
The frequencies and damping time of the quasi- normal oscillations called " quasi- normal 
modes" (QNMs) are determined only by the black hole's parameters and are independent 
of the initial perturbations. It should be possible to infer the black hole parameters solely 
from the QNM frequencies. QNMs carry information of black holes and neutron stars, and 
thus are also of great importance to gravitational- wave astronomy These oscillations are 
mainly produced during the formation phase of compact stellar objects and can be strong 
enough to be detected by several large gravitational wave detectors under construction. Re- 
cently, QNMs of particles with different spins in black hole spacetimes have also received 
much attention 

Numerical relativity has provided evidence that the QNMs dominate the gravitational- 
wave signal associated with many processes involving dynamical black holes (such as the 
formation of black holes in gravitational collapse or binary merger). Since the QNMs encode 
information concerning the parameters of the black hole one may hope that the gravitational- 
wave detectors that are now coming into operation will be able to use these signals to inves- 
tigate the black-hole population of the universe. Even though most studies of QNMs have 
been motivated by their potential astrophysical relevance, there are several other reasons 
why one might be interested in understanding the spectrum of oscillations of a black hole. 
In particular, the modes have played a role in discussions of black-hole stability |3(] . Recent 
investigations of large extra dimensions, has led to the somewhat strikingprediction that 
mini black holes may be observed at particle accelerators such as the LHC|4]]. One problem 
is that in order to suppress a rapid proton decay we need to physically split the quarks 
and leptons. Such models are generically called split fermion models[5]. Most methods in 
evaluating the QNMs are numerical. Recently, the Dirac field QNMs were evaluated for a 
Schwarzschild black hole [6]. A powerful WKB scheme was devised by Schutz and Willj3], 
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and it was also extended to higher orders [3]. In particular, the interaction of fermions with 
a gravitational field is an important topic which was firstly studied by BriU and Whee.erS 
in the 1950s, and later by many other authorsflOj. It is well known that the vector field 
equation in a Schwarzschild spacetime can be separated and the time and angular equation 
integrated. One is then left with the separated radial equation whose integration is not im- 
mediate. The remaining Schrodinger-like equation has received great attention and provides 
QNMs. 

On the other hand, the neutrino does not respond directly to electric or magnetic fields. 
Therefore, if one wishes to influence its orbit by forces subject to simple analysis, one has 
to make use of gravitational fields. In other words, one has to consider the physics of a 
neutrino in a curved metric. In this work we revisit this topic using some ideas of the Fock- 
Ivanenko formalism and SUSY to study the radiation of energy produced by a neutrino, 
which is affected by a gravitational field produced by a Schwarzschild BH. To do so it we 
shall study approximated analytical solutions for the Dirac equation for a massless 1/2-spin- 
fermion around a neutral and non-rotating Schwarzschild BH. This topic has been recently 
studied using numerical methods, in a d + 1 spacetime, where the d — 2 extra dimensions 
are compact jllj. In this paper we expand the formalism there developed, making use of the 
supersymmetrical (SUSY) properties of the Schrodinger equation. 

The paper is organized as follows: in Sect. II we obtain the Dirac equation for a 1/2-spin 
massless fermion without charge (neutrino) in a static and spherically symmetric spacetime, 
using the Fock-Ivanenko formalism. In Sect. Ill we obtain the radial equation for the spinors 
on such a metric. The spinors obey Schrodinger-like stationary equations. In particular, in 
Sect. IV we study these equations for a Schwarzschild BH. In Sect. V we study the SUSY 
properties of the effective radial potentials of the Schrodinger-like equations. In Sect. VI 
we obtain approximated solutions of the quasi-normal frequencies and their eigenstates and 
SUSY partners for the spinors. In both asymptotic cases, he neutrinos are close, and very 
distant to the BH. The quantization of quasi-normal modes is studied in the appendix. In 
Sect. VII we have included some comments about the treatment of d > 4- neutrinos where 
the additional dimensions are compact. Finally, in Sect. VIII we develop the final comments. 
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II. DIRAC EQUATION FOR A NEUTRINOS IN A STATIC AND SPHERICALLY 
SYMMETRIC SPACETIME 



We consider the metric that describes a static and spherically symmetric spacetime 

ds 2 = f{r)dt 2 - f- l (r)dr 2 + r 2 dn\ (1) 

where Q is the solid angle, such that dfl 2 = sin 2 (#) dO 2 + cos 2 (#) d(p 2 . In this paper we shall 
use natural units: H — c — 1. The Dirac equation for a free massless fermion \1/ with spin 
1/2, on the metric ([I]), is 

f%# = 0, (2) 

where V M denotes the /ith component of the covariant derivative and 7 M are the Dirac matrix 
components on the metric ([1]). Using conformal transformations, one obtains 

9w -> 9nv = Q 2 9iiu, (3) 
f -> * = 0^ 3 / 2 ^, (4) 

7"V M # e 5 / 2 7 M v^, (5) 



where we shall take O = 1/r, so that ^ = r 3 / 2 \I/. The resulting metric is 

ds 2 = ^dt 2 - -±-dr 2 - dVL 2 , (6) 



on which 



7„V"* = 0. (7) 



The spinors x {fii 0) on the 2-sphere obey the equation 1 12 j 



%V> X ™=i I ]x, <tW =i(«tJx, (n) , (8) 



where 



«tt= I ^ I • ( 9 ) 



1 Notation: in what follows we shall use the compact notation denoting by t all that correspond with the 
spinor "up", and with \. all that correspond with the spinor "down". For instance, the scalar = I + 1 
corresponds to x an d Kj, = —2 corresponds to x ■ 
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Notice that (/Cf^) can take integer values, but the zero is excluded: (k^) 7^ 0. If we expand 
the function as 

• -E^xf+W), do) 



the equations of motion for 0^ and 0^ as a result are found to be[13] 



{T^o + T'Vi + i (/ + 1)7 5 } 
{7°V + 7^1-^7^} 



51 ,(t) 

5/1 (40 



0, J = 0,1 



0, 1 = 1,2,..., 



(11) 
(12) 



where we use the spatiotemporal coordinates (t, r, 9, 0) and 7 s = 7°7 1 7 2 7 3 is the interaction. 
With the aim to solve the equations fill I) and (I12p . we shall make use of the Fock-Ivanenko 
formalism 14| . To represent the system we choose the following Dirac matrices: 

<yt = -t—L= a 3 , Y = r^/J{r)a 2 . 

v/( r ) 

where the o % are the Pauli matrices 



a 




a' 




a" 




so that the spin connections, V ^ = jj^^-^rg^ 2 , are given by 



T> = a l 



l (r 2 \ d ( f(r) 



4 / dr \ r 



0. 



(13) 



(14) 



(15) 



III. EQUATIONS OF MOTION FOR 



(t) 



A. Spinor <fi) 

The equation of motion for 0^ (the equation for 0^ can be worked out in the same 
way) is 



d 

77" + 



T d I Vf( r ) 



dr 2y/f(r) dr \ r 



(t) 



i(t) 



V 7 ^ 



dt 



2 The reader can see, for example, page 113 - chapter 6 - of the book 

3 Notice that the sum does not run over /z, but it does over v. 
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We propose solutions of the form 



z (t) (r,t) 



a/70 



-1/2 



n v I 



(17) 



where is some constant of integration. With this transformation, the Dirac equation 



can be written as two coupled first-order equations: 



f(r) 

f(r) 



(Z + l)*S f) (r) = E*>*™{r] 
(Z + l)*S f) (r) = -^PM t} ( 



(t)xT/W ( 



fir 



+ 



dr r 

Now we define the tortoise coordinate u and the function W(r), as 



ar aw r 



so that the equations ( fT8"j) and ( TIT?]) can be written, respectively, as 



A_ W (t)( M ) 
aw 

A + wr(t)( u ) 
du 



These equations can be separated as 

d 2 



du 2 
du 2 



+v 1 %)J M t} («) = W } ] 



(18) 
(19) 

(20) 

(21) 
(22) 

(23) 
(24) 



such that the supersymmetric (SUSY) potentials and V^ 1 ' are given by the expression 



(25) 



The effective potentials in (T25]) are supersymmetric. Hence, the functions and ^ have 
the same spectrum for quasi-normal modes. 



B. Spinor ffl 

In the same manner one can obtain the Schrodinger-like equations once we propose the 
solutions 



<^(r,t) 



-1/2 



-iE, 



*£»(r) 



(26) 



such that the superpotential W™(r) is 

W^\r) = - l -^m, (27) 
and the Schrodinger-like equations for ^f^ ^ are 




with the potentials V-y\ (u) given by 

V$(n) = r^+[wW(n)] 2 . (30) 

The equations (I23H24P and (I28H29P give us the radial information about the wave functions of 
the spinors (f)\ (r,t). In the following sections we shall consider an analytical treatment to 
solve the pairs of Schrodinger-like equations ( 1231124")) and (12811291) . for a Schwarzschild black 
hole. 



IV. NEUTRINOS IN A SCHWARZSCHILD BLACK-HOLE 



We consider the particular case where the function f(r) in the metric is given by 



/(r) = l--, 
r 



(31) 



such that a = 2GM is the Schwarzschild radius. In this case the tortoise coordinate u(r) is 
given by 



- 1 



or, 



u(r) = r + a In 
From the second equation in ( 1321) . we obtain 

[1 " /(«)] 



u(r) 



1 - f(r) 



a In 



f(r) 



Li -/(»•) 



(32) 



(33) 



1 + LW ( e («- a )/ a )' 

where LW(x) is the Lambert function, such that LW(x) e LW ^ = x and LW(0) = 0. Now 
we are aimed to obtain the potentials V^\u) as a function of u. From eqs. ( 1201) . the second 
equation of fl32|) and fl33|) . we obtain the superpotentials 

(n n ) [W[eM/-]]^ 



iy (n) ( 



u 



[1 + LW [ e («^)/«]] 3/2 ' 



(34) 



The exact potentials, written function of u, are 



du 



LW [e^-^l 11 ] 



[1 + LW [e^ u - a )/ a }\ 



2 a 2 [i + LW [ e («-«)/a]]7/2 V 



(35) 



In the limit case where the particle is close to the Schwarzschild radius: u — > — oo, so that 
e («-o)/o _j, q j n cage e ff ec tive potentials can be approximated by 

>'-' . . (,. \ . . .. . 

(36) 



V in) (v) ~ fa-aVa , ( K tJ (ii-a)/(2o) 

^1,2 W " .2 e ± 2a 2 



In order to make a dimensionless description of the problem we can make the following 
transformation: 2v = —(it — a)/a, and the Schrodinger-like equations ( |23i) and f l24l) . for 
(u) and ^3(1;), close the Schwarzschild radius, become 

2 



d 2 



dv 2 



E. 



E 




0, (37) 



where E$ n) 



E%' n) and 



1, j = 2 a = n + 1, (3 = n, for 

2, j = 1 a = n, /3 = n + 1, for 



(38) 
(39) 



Notice that for large r one obtains u — > — oo, so that V^ lf — 0- I* 1 other words, the 
approximated potential (1361) becomes null for large distances and takes the same asymptotic 
value as the exact potential fl35|) . Therefore, the asymptotic solutions of the eq. (1371) will be 
a good approximation for r — > oo. However, the approximated potential (1361) is not a good 
approximation for intermediate distances. 



V. SUPERSYMMETRY OF HAMILTONIANS 

We consider the Hamiltonians and 



^ = -^ + [^']'W-^, (40) 



H™ = -^ + [^>]>) + ^, (4!) 
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such that l4 n) 



Hi 



(U) _ dW^(v) 



dv 



. In our case W^(v) ~ — 2(/^)e vi , so that we obtain 



— = 2W /( -^- ) (f ). Following Mielnik 15], we shall define the differential operators: 



r(U) 



dv 
dv 



(42) 
(43) 



then, we see that Hamiltonians f|4*0|) and (j4"Tj) can be rewritten in terms of b(^) and bj^: 

b(u)b ( n) = hwHn) + Hn)Mn)] = hnWn) + ™^\v) = H<™ + 2W^\v) = Hp- 

(44) 



(til 

Furthermore, H; b 



Therefore, if ¥? l \aJ, v) is and eigenvector of FF- with eigenvalue 



m> : b (n) b (t4-) b (n) 



(tio 



b 



which means that 



2 -W>) 



b (T4) W j 



(45) 



(n) 



= b*^^ 1 w (/3, Z, u) will be an eigenvector of Hj-'^ with eigenvalue 
In other words the potentials and V„- possess the same spectra of quasi-normal fre- 
quencies (but with different eigenstates), because they are supersymmetric partners derived 
from the same superpotentials W^(v). Notice that in cases where En is not real, the 
Hamiltonians will not be hermitian. 



,(t4-), 



(in) 



hence 

! 1 

n 



p(U) 

J-"l 



VI. APPROXIMATED SOLUTIONS 

In order to solve the equations (|37|) we must consider that such equations are only valid 
close to the Schwarzschild horizon: r > a. In this limit, v — > oo. We shall consider separately 
the cases for the potentials and V 2 ; , for the spinors <f>\ ; . 

ft) 

A. Spinor ' 

As a first case we consider the radial equations for the spinor <p\ . We shall consider, 
separately, the attractive and repulsive potentials in the equation (1251) . 

4 Actually, there exist a family of solutions W^(v) = g^tJ^Ei^-^ for Ei ( az ) _ 

z a-i p(x _ flj z ^ w hich give us SUSY potentials v}J 2 l) (v) = 4(n n ) 2 e~ 2v ± 2(K n )e~ 1 \ but in this paper we 
shall restrict our study to the particular solution with C = 0. 
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1. Attractive potential for ^ 



(t) 



The first case consists in 



V®\l, v) = 4(1 + lfe-' 2v - 2(1 + 1) e~ v 



(46) 



which is similar (but not exactly equal) to a Morse potential 5 . The general solution for A3 Tj) 
is given by 

dv e v 



^\nM~e^ 2 [l v A\z(v)]-Z V2 [z(v)]] 



B n l + C n l 



[Z^ -2^, [«(«)]] 

where the Z^[z(u)] and X u , 2 [z(v)] are the modified Bessel functions of the first kind, with 



(47) 



^1,2 



2iaE^^~ 



z ( v ) = 2(1 + 1) e~ v . 
As can be demonstrated [see appendix (jA"]) ] the eigenvalues of energy are 6 : 



(48) 
(49) 

(50) 



and the eigenvalues for the angular moment are I > 0. In order to avoid a possible divergence 
in ( 1471) . we choose C n \ = 0. A particular manifestation of Hamiltonians which are SUSY, is 
related to the null eigenvalue of energy: E^l = 0. In this case n = and the solution is 



*$\0,l,v) ~ S w e-" /a [X_i /2 [«(i;)] - X 1/2 [z(«)]] oc e 



-2(J+l)e-° 



(51) 



where Bqi is a constant to be determined by normalization. 

The \l/ 2 ( n 5 ^ f)-partner states can be obtained when we apply the operator of fej: to 

¥?\n + l,l,v) = |-A + 2 (Z + l)e-"| ^(n,/^). (52) 

The correspondence (1521) is a manifestation of the SUSY character of the Hamiltonians 
and H-P. 



It is interesting notice that the original bound state Morse potential has the form VM orse (v) 



S 2 [e- 



2 c 



and corresponds to an Hermitian exactly solvable potential with real eigenvalues of en- 



ergy. However, as was demonstrated in (19) . it is possible to extend the usual theory for exactly solvable 
(ES) 2C§ and quasi-exactly solvable 2l| (QES) potentials to accommodate QNMs solutions. 
6 In general the energy can take the values — =pm/(2a), but we choose only the eigenvalues with sign 
minus because they are those that correspond to decaying modes. 
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2. Repulsive potential for 



(t) 



Now we consider the potential 

vP(v) = 4(Z + l) 2 e~ 2v + 2(1 + 1) e~ v . 
In this case the general solution is given by 



(53) 



¥?\n + l,l,v) ~ B 



ni 



( i + l) e -«/2 + ( n+ _) e «/2 



X w [^)]+ e -/ 2 (Z + l) 



+ a 



(n + l)4 n+1 16~ (n+1) 



4-(»+i) 



e «/ 2 + e -«/2 (/ + i) 4 («+i)i 6 -(«+i) 



+ e -«/ 2 4( w+1 >16- (,l+1) (Z + l)X^[z{v)]} 
If dve v 



x 



where 



A J { [ Z + i + ( n + i) e „] j m [^)] + (/ + 1)Xm2 [^)] } 

)Ui = n + l/2, /i 2 = ^ + 3/2, 
z( v ) = 2(Z + 1) e^. 



2 ■ 



^1 [«(«)] 



(54) 



(55) 
(56) 



In order to the SUSY expression for the partners (1521) to be fulfilled, we shall require that 
_B n ; = 0, so that the resulting nonzero constants B n i in (j!7|) and (7 n z in (154]) should be deter- 
mined by normalization. As can be demonstrated [see appendix (TA~]) ]. the energy eigenvalues 
are 

^-^-^ » = ^-- (-) 

Notice that M 2,t) = M 1,t) [see eq. (Ej 



B. Spinor 



For completeness, we consider the radial equations for the spinor such that the 
attractive and repulsive potentials are given by ( 130]) . 



1. Attractive potential for 



Now we consider the attractive potential related to = —I. In this case the potential 



is 



V} i] (v) =4l 2 e-' 2v -2le~ v , 



(58) 
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so that the general solution of (157)1 is 

*?(n,l,v) ~ e"^ 2 {a ni [X n _ 1/2 [-2Ze~*] + X n+1/2 [-21 e^]] (59) 
+ /?n* [/C„-i/2 [-21 e- v ] - /C n+1/2 [-2/ e~ v ] } } , (60) 

such that the functions /C nT i/ 2 are the modified Bessel functions of second kind, and the 
quasi-normal frequencies are 

E^ = -i^, n = 0,l,.... (61) 

The particular case with n = give us the bottom energy eigenvalue E^j = 0, which 
corresponds to the eigenfunction 

*?(0,l,v) = a 0l e- 2le - v , (62) 
where we have put (3 n i = 0, for the eigenfunction to be finite along all the domain v. 

2. Repulsive potential for 0j 

Finally, for completeness, we consider the repulsive potential associated to tv>^ = —I. 
The potential has the form 

V 2 il \v) = Al 2 e~ 2v + 2le~ v , (63) 
and the general solution in this case is 

^\n+l,l,v) ~ e' v ' 2 {a nl [X n _ 1/2 [21 e^] +X n+1/2 [2le~ v ]] 

+ 0ni [lCn-i/2 [21 e- v ] - K n+1/2 [21 e' v ] } } . (64) 

Finally, the quasi-normal frequencies are 

E^ = -i^ n = l,2,.... (65) 

In order to fulfill the SUSY expression for the partners 

^\l,l, v ) = !-±-2leA *?(0,l,v), (66) 



dv 

we shall require that a i = and more generally, that a n i = in (|6 
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C. Asymptotic solutions for large distances: r — > oo 



To complete our study we shall study the weak gravitational field case for large dis- 
tances to the BH. In this case r — > oo, so that ^ u - a )l a _ ). e"/ a . The superpotential can be 
approximated by 

w {n) (u)\ i m 

W \ Vr ^ D - a LW { e u/ a y 
so that the SUSY potentials will be 

V$\u)\ ' , IMT11. (68) 



— > 0, in agreement with 



For very large distances these potentials tend to 0: vffi (u) 
the approximated potential ( 136]) . Hence, the approximated Schrodinger-like equations for 
u — t- oo can be written as 



d 



2 



^-^(u)^[E^] 2 ¥^(u), (69) 



du 



(1*40 

which has a general solution that can be written as a linear combination of e ±lEn u . After 
taking into account the normalization conditions, and the signature of E n , we obtain that 
for large distances the outgoing solutions are 

~ e - nu/(2a) . (70) 

This solution agrees perfectly with the whole obtained in 22]. 



D. Quasi-normal frequencies and Hawking temperature 

The Hawking temperature for the Schwarzschild BH is 

1 



H 



Aira 

so that the quasi-normal frequencies can be written as 



(71) 



= Eg>™ = -2mnT H . (72) 



This result is exactly whole obtained recently in 111], but using the WKB method (see 
also[16|). 



13 



VII. NEUTRINOS IN A d > 4-DIMENSIONAL SCHWARZSCHILD BH 



The study of neutrinos which are close to multidimensional Schwarzschild BH is an in- 
teresting issue. In this case the extended d > 4-dimensional Schwarzschild metric is given 
by[23| 

ds 2 = f [d) (r)dt 2 - f- d ]{r) dr 2 + r 2 dtf d _ 2 , (73) 

where d is the dimension of the spacetime, dVt 2 d _ 2 denotes a metric of a (d — 2)-dimensional 
sphere and 

/ M « = l-0 ( "' (74) 
In this case the scalars analogous to (J5]) are 



K 



(d) 

n 



(75) 



where I > -(d - 2)/2 + 1 and / > -(d - 2)/2 + 2 for n ( f and k^> , respectively. Note that 
d only can take integer-even values in order for k$ to be a nonzero integer. The tortoise 
coordinate can be written as 

(3-d) 



u(r) = a 



I I ,/- ! 

21 < L '3-rf'rf-3 ; 



i-/(«o(0. 



l/(d-3) 



1 



(76) 



where 2 ^i[/i, ^; 7; = En 
fact that 



OO Un^nX 

=0 7„n! 



is the hypergeometric function and we have used the 



(a/r) 



1 



l/(d-3) 



In this case one can approximate the superpotential for very large extra dimensions, by 



(77) 



W {n{d) \r) 



d—^oo 



l/(3-d) 





/U \l/(3-«0" 


(3-d) - 


1 - 


. 1+ (i-0 











1 1/2 



(78) 



The detailed analysis in the SUSY framework is very complicated and goes beyond the scope 
of this work, but finally one can find that the eigenvalues of the energies are 22 ] 

'd-3 s 



E (n(d)) _ _ t 



2a 



n = — 2tt iTn n, 



(79) 



so that the d-dimensional Hawking temperature increases linearly with the number of extra 
dimensions: T# = (d — 3)/(47ra). 
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VIII. FINAL COMMENTS 



In this paper we have used the Fock-Ivanenko formalism for the Dirac equation in curved 
spacetime to write the Dirac equation for a massless 1/2-spin-fermion around a Schwarzschild 
BH. We have calculated approximated analytical solutions for the spinors. The radial eigen- 
functions of the spinors <p\ are described by Schrodinger -like equations for the radial 
eigenfunctions of the spinors, which can be rewritten in terms of the redefined tortoise co- 
ordinate u(r): v = —{u — a) / a. We have proven the SUSY character of the potentials 
and V 2 f , corresponding to the Hamiltonians and H 2 , which are not hermitian. 
We obtained approximated analytical solutions for eigenvalues of energy (or quasi-normal 
frequencies). The interesting result here obtained is that all the asymptotic states (and 
their SUSY partners) have purely imaginary frequencies as eigenvalues of energy, which can 
be expressed in terms of the Hawking temperature T#: E n = —IninTn. Furthermore, 
the isolated bottom state has a real null eigenvalue of energy, in agreement with what one 
expects for a SUSY Hamiltonian 24J. Therefore, all the asymptotic modes decay, except the 
whole with zero energy. The treatment with d > 4- compact extra dimensions deserves a 
more detailed study and goes beyond the scope of this work. However, some comments were 
included in Sect. VII. 
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Appendix A: Quantization of quasi-normal frequencies: states 

We consider the equation ( |37|) for the massless Dirac modes ^2 ■ The 

case with = 

{£ + 1) accounts for the ^jj, (with £ = 0, 1,2, ...), while k± = —£ corresponds to the ^/[^ 
ones (of course with £ = 1, 2, ...). 

Although two elections are possible for the potentials (the second exponential term in 
the potential can take signs + or — ), the Schrodinger equation itself seems to come from a 
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Morse-like potential problem. In symbols: 

- + {^n) 2 e~ 2v + 2(K U )e-v - 4a'[E£>™]*}9™ = 0, (Al) 

+ {4(^u) 2 e- 2 ^ - 2(Kt,)e- - Aa 2 [E^} 2 }^ = 0. (A2) 

Then, we take the way followed by Morse to solve our problem First, we make a 
change of variables: 

- v d d d 2 d 2 d 2 

^ ' dv ^ dy' dv 2 ^ dy ^ dy 2 

The resultant equations are 

d 2^(n) ld9 (u) [e (i,u) ]2 2( ) 

-ry- + ^- + " 4(«tJ }*i = 0, (A4) 

f + + 4(Kn)2} ,«u» = o, (A5) 

dy 2 y dy y 2 y 

where the parameters 4a 2 [En'^] 2 have been replaced by [£n ] 2 . Finally, a clever transfor- 
mation of functions bring us closer to a problem with hypergeometric differential equations. 
Let the set of functions be F^>^(y) such that 

^\y) = e- a y(2ay^Fi U \y), (A6) 

^ i \y) = e~ a y{2ay)^F^\y). (A7) 
At the same time we define the variable z as 

z = 2ay, A = 2 «A jL = Aa 2 ^-. (A8) 
dy dz dy 2 dz 2. 

Putting it all together, we arrive at the final expression, which entails particular cases of the 
general confluent hypergeometric differential equation [131] 

d 2 F[ n) „ a M n) 
1 + (l + p-z^ 1 



dz 2 dz 

cM.. \ i ../fl i 1\ _.2 /I/.. \2 "I 



4[^F+^ 2(« t J + q(£ + l) a 2 -4(^) 2 
^ 4i 2a + 4^ ^ " U ' 



d 2 F^ . a M n) 



z^r- + (1 + P-z 



dz 2 1 dz 

rt/.. \ _Y o i 1\ _.2 /I/.. \2 "I 



4[^] 2 + /3 2 , 2(« n )-a(£ + l) a 2 -^) 2 ,,, 
' ^ 4^ + 2a + Aa 2 Z <* 2 
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Under certain conditions (which we shall study later), these equations have the form of a 
confluent hypergeometric differential equation 

,fMl + {1 + p- z) ^M + nF ^ t) = . (A9) 

The general solution can be expressed in terms of the Laguerre polynomials L(n, ft, z) 

V{ft) L(n + ft,-ft,z) , T(-ft) L(n,ft,z) 



A n/3 L(n,ft,z)+B n £ 



, (A10) 



r(— n) bin(n, ft) zt 3 r(— n — ft) bin(n + ft, n) 

where bin(n, ft) — n\j [ft\(n — ft)\]. Whatever the rest of the arguments, if n becomes positive 
integer, the solution always will be able to be written as a polynomial form (finite number 
of terms). This treatment of the solutions is particularly useful for the states; and we have 
the ground states (n = 0) using the fact that F 0/3 = const, in the attractive potentials. 
There are four significant cases of interest, which we shall study separately. 



a. K t = I + 1 and F x (t) 
The first case is 

4[ei U) ] 2 + /3 2 = 0, (All) 
2(K t )+a(ft + l) 



(A12) 
2a 

a 2 -4(/« t ) 2 = 0. (A13) 

A decaying exponential factor in f ]A6j) forces us to set a = 2(1 + 1); then we obtain ft = 
—2(n + 1). Finally £n ^ = = —i(n + 1). The minus election in the sign of £ is because 



of the energy 



E M) = = -i^H±3. (A14) 
2a 2a K J 



which must have a negative imaginary part to recover the outgoing quasi-normal modes. 
Notice that n must be a positive integer, and due to the fact ft = —2{n + 1) < 0, we obtain 
n = 0,1,2,3,... 
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K t =£ + 1 and F 2 (t) 



In this case 



4[e t) ] 2 + /? 2 = 0, (A15) 
feL^tl! = »', (A16) 
a 2 - 4(/t t ) 2 = 0. (A17) 



(2,t) _ .0 



?7T 



Using similar arguments, but with <\A7\i . we obtain a = 2(£ + 1), (3 = —2n' and £ 
—in'. Again n' must be a positive integer. 

E&V = ^_ = (A18) 
2a 2a K ' 

Using the fact that f3 = —2n' < 0, we obtain n' = 1,2,3, so that n = n' — 1, and the 
energies ( 1A18I) can be written as 



n(l,t) = ? - ( n + 1 ) 

where n = 0, 1 



= E M) = (A19) 



c. K; = -| and 



Only the main results are of interest for us; a = 2£, (3 = —2m and £m = i§ = —im, 



2 



with m a positive integer. The energy in this case is 



£(U) = ^_ = (A20) 
2a 2a V ; 



where m = 1, 2, ... 



d. k± = -£ and F 2 a) 

In this case a = 2£, (3 = —2{wl + 1) and £m = if = —i(m' + 1), with m' a positive 
integer. The energy in this case is 
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A similar analysis to those made in (1A0 aj) and (lAObl) , gives us m! = m — 1 and finally we 
find that the energies (or quasi-normal frequencies) have a unique result (we relabel m — > n) 

£(u) = E (U) = _ i JL j (A22) 

where n — 1,2, .... 
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